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We present the results of a search for the rare flavor-changing neutral-current decays B — > n£ + £~ 
(it — tt ± ,tt° and 1 = e,fj) and B° — > -n£ + l~ using a sample of e + e~ — s> T(4S) — ¥ BB decays corre- 
sponding to 428 fb _1 of integrated luminosity collected by the BABAR detector. No significant signal 
is observed, and we set an upper limit on the isospin and lepton-flavor averaged branching fraction of 
B(B tt£ + £~)< 7.0xl0~ 8 and a lepton-flavor averaged upper limit of B(B° ->■ r)£ + l~)< 9.2xl0~ 8 , 
both at the 90% confidence level. We also report 90% confidence level branching fraction upper 
limits for the individual modes B + — ¥ ir + e + e~, B° — > 7r°e + e~, B + — > 7r + /i + /i~, B a — > 7r°/i + /i _ , 
B° -> ne + e~, and B° — > rjfi + fi~ . 



PACS numbers: 13. 20. He, 13.60.Hb 
I. INTRODUCTION 



In the standard model (SM), the decays B — > ir£ + £~ 
(tt = 7T ± ,7r° and I = e,/i) and B° —> r/£ + £~ proceed 
through the quark-level flavor-changing neutral current 
(FCNC) process b -> d£+£~. Since all FCNC processes 
are forbidden at tree level in the SM, the lowest order dia- 
grams representing these transitions must involve loops. 
For b — > d£ + £~, these are the electroweak semileptonic 
penguin diagrams (Fig. 1(a)) and the W + W~ box dia- 
grams (Fig. 1(b)). The b —> d£ + £~ transition is similar 
to b —> s£ + £~, but its rate is suppressed by the ratio 
|Vjrf/T4s | 2 ~ 0.04 where V t d and V ts are elements of the 
Cabibbo-Kobayashi-Maskawa quark mixing matrix 
The predicted branching fractions for the B + — » ir + £ + £~ 
decay modes lie in the range of (1.4-3.3) x 10~ 8 , when 
the dilepton mass regions near the J/tjj and ijj(2S) are ex- 
cluded in order to remove decays that proceed through 
the intermediate charmonium resonances. The largest 
source of uncertainty in these predictions arises from 
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knowledge of the B — >• n form factors These 
branching fractions imply that 5-15 events occur for each 
B -> tt£ + £~ decay channel in the BaBAR data set (471 
million BB pairs). The predicted B° — > r/£ + £~ branching 
fractions lie in the range (2.5-3.7) x 10~ 8 where again 
the dominant uncertainty is due to lack of knowledge of 
the B° — > r\ form factors 0]. 

Many extensions of the SM predict the existence of 
new, heavy particles which couple to the SM fcrmions and 
bosons. The b — > d£ + £~ and b — s- s£ + £~ decays provide 
a promising avenue in which to search for New Physics 
(NP). Amplitudes from these NP contributions can inter- 
fere with those from the SM, altering physical observables 
(e.g., decay rates, CP, isospin, and forward-backward 
asymmetries) from the SM predictions 0, 0, @~G3 ■ Mea- 
surements in the tt£ + £~ and rj£ + £~ systems complement 
and provide independent searches of NP from those in 
the K.(*H + £~ channels [llT - |l~7l |. as physics beyond the 
SM may have non-trivial flavor couplings [rsj . The mea- 
surement of observables as a function of the square of the 
invariant dilepton mass q 2 = m 2 e for exclusive b — > d£ + £~ 
decay modes allows for more thorough tests of SM pre- 
dictions and deeper probes for NP but is currently not 
possible due to the size of the data set. 

Only one b — > d£ + £~ decay has been observed to date, 
with LHCb measuring the B + — > 7T + (i + [i~ branc hing 
fraction to be (2.4±0.6±0.2) x 10" 8 [H . Both BABAR [if 
and Belle [2l[ have performed searches for B — > n£ + £~ 
decays, but have observed no significant signal. For the 
tt£ + £~ modes, the smallest upper limits from the B fac- 
tories lie within an order of magnitude of the SM pre- 



6 




FIG. 1: Lowest order Feynman diagrams describing the quark level b — > d£ + £ transition in B meson decay: (a) electroweak 
penguin diagrams and (b) W + W~ box diagrams. 



dictions [3|-|5| and are beginning to exclude portions of 
the NP parameter space. No previous searches for B° —> 
r\£ + £~ have been reported. Observation of b — > d£ + £~ de- 
cays at the B factories is currently limited by the size of 
the available data sets. Additionally, for B + — > it + £ + £~ , 
background from B + — > K + £ + £~ decays where the kaon 
is misidcntificd as a pion must be treated carefully as 
K + £ + £~ can appear very signal-like and occurs at a rate 
approximately 25 times the expected B + — > n + £ + £~ra,tc. 

In this article we report on our study of the B — > 
ir£ + £~ and B° -> r]£ + £~ decays using the full BABAR data 
set. presenting branching fraction upper limits for the 
modes B + 7T + e+e~, B° -> n°e + e~, B + -4 n + //+ , 
B° 7t°fi, + fi~, B° -> r/e+e", and B° -> n~ . Charge 
conjugation is implied throughout unless specified other- 
wise. We also present upper limits for the lepton-flavor 
averaged modes B+ -> n+£+£~, B° -> tt ^", and 
B° — > r]£ + £~, where we constrain the e + e~ and 
branching fractions to be equal; for the isospin averaged 
modes B — > iie + e~ and B — > ir/i + [i~, where the B + — > 
ir + £ + £~ decay rate is constrained to be twice the B° —> 
tt°£ + £~ decay rate; and for the isospin and lepton-flavor 
averaged mode B — > tt(' + £~. For the lepton-flavor av- 
eraged measurements, we neglect differences in available 
phase space due to the difference between the electron 
and muon masses. The branching fractions are based 
on signal yields that are extracted through an unbinned, 
extended maximum likelihood fit to two kinematic vari- 
ables. All selection criteria are determined before the fit 
was performed on data, i.e., the analysis is performed 
"blind". 



II. BABAR DETECTOR, SIMULATION, AND 
DATA SETS 

The results of this analysis are based upon a sample 
of e + e~ — > T(4S) — > BB interactions provided by the 
PEP-II asymmetric-energy storage rings and collected by 
the BABAR detector located at SLAC National Accelera- 
tor Laboratory. The BABAR data sample corresponds to 
an integrated luminosity of 428 fb -1 containing 471 mil- 



lion BB decays. This is the full data set collected at the 
T(4S) resonance. A detailed description of the BABAR de- 
tector can be found elsewhere [22| . Charged particle mo- 
menta are measured with a five-layer, double-sided silicon 
vertex tracker and a 40-layer drift chamber operated in 
proportional mode. These two tracking systems are im- 
mersed in the 1.5 T magnetic field of a superconducting 
solenoid. A ring-imaging Cherenkov detector with fused 
silica radiators, aided by ionization loss dE/dx measure- 
ments from the tracking system, provides identification 
of charged particles. Electromagnetic showers from elec- 
trons and photons are detected with an electromagnetic 
calorimeter (EMC) constructed from a finely segmented 
array of thallium-doped Csl scintillating crystals. The 
steel flux return of the solenoid (IFR) was initially in- 
strumented with resistive plate chambers (RPCs) and 
functions primarily to identify muons. For the later data 
taking periods, the RPCs of the IFR were replaced with 
limited streamer tubes and brass to increase absorption. 

The BABAR Monte Carlo (MC) simulation utilizes 
the Geant4 package [23| for detector simulation, and 
EvtGen [HI and Jetset7.4 [H] for BB and e + e~ qq 
(q = u,d,s,c) decays, respectively. The BB and contin- 
uum (e + e _ — > qq, q = u,d,s,c) MC samples correspond 
to an effective luminosity of about ten times the data 
sample collected at the Y(AS) resonance. Simulated B 
— > n£ + £~~ signal decay samples are generated according 
to the form- factor model of Ref. [2611 . with the Wilson 
coefficients taken from Refs. |27H29|. and the decay am- 
plitudes calculated in Ref. [10| . The B° — > rj£ + £~ signal 
MC samples utilize the same kinematics, Wilson coeffi- 
cients, and form-factor model as the ir£ + £~ modes. The 
effects of the choice of form-factor model and the val- 
ues of the Wilson coefficients are considered as sources 
of systematic uncertainty in the signal efficiency. We 
also make use of simulated B -> K^£ + £~, B -> J/ipX, 
and B — > ip{2S)X samples. Signal efficiencies, as well 
as parameters of the fit model, are determined from sig- 
nal and K^*H + £- MC data sets. The B — » J/ijjX and 
B %[)(2S)X MC samples allow us to study background 
from these decays and also serve as the data sets from 
which we fix the parameters of the fit model used in the 
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B J/i/jw/r] fit validation, as described later. 



III. EVENT RECONSTRUCTION AND 
CANDIDATE SELECTION 

Event reconstruction begins by building dilepton can- 
didates from two leptons (e + e~ or /u + /i~). Leptons are 
selected as charged tracks with momenta in the labo- 
ratory reference (Lab) frame greater than 300 MeV/c. 
Loose particle identification (PID) requirements are 
placed upon the two leptons. More stringent PID require- 
ments are applied later, and the optimization of these se- 
lection criteria is discussed in Section [V] The lepton pair 
is fit to a common vertex to form a dilepton candidate 
with the requirement that mu < 5.0GeV/c 2 [30fl . We 
also place a loose constraint on the x 2 probability of the 
vertex fit by requiring it to be greater than 10 -10 . For 
electrons, we apply an algorithm which associates pho- 
tons with electron candidates in an attempt to recover 
energy lost through bremsstrahlung, allowing at most one 
photon to be associated with each electron. The photon 
trajectory is required to lie within a small cone of open- 
ing angle 0.035 rad about the initial momentum vector 
of the electron, and the photon energy in the Lab frame 
must be greater than 30 McV. Additionally, we suppress 
background from photon conversions by requiring that 
the invariant mass of the electron (or positron) paired 
with any other oppositely charged track in the event be 
greater than 30 McV/ c 2 . 

Charged pion candidates arc charged tracks passing 
pion PID requirements which retain approximately 90- 
95% of charged pions and only 2-5% of charged kaons. We 
reconstruct ir° candidates from two photons with invari- 
ant diphoton mass m 77 lying in the range 115 < m 77 < 
150MeV/c 2 . A minimum value of 50MeV is required for 
the Lab energy of each photon. Photons are detected as 
EMC clusters not associated with a charged track. The 
clusters are also required to have a lateral shower profile 
consistent with originating from a photon. We recon- 
struct rj as 77 — > 77 (?y 77 ) and 77 — > 7r + 7r~7r° (773^), which 
constitute 39.3% and 22.7% of the 77 branching fraction, 
respectively. As in the case of the 7r° , we require the ?y 77 
photon daughters to have energy greater than 50 MeV in 
the Lab frame. Additionally, the photon energy asymme- 
try A 1 — \E\ tl — E2 n \/ {E\ n + -£2,7) must be less than 
0.8, where E\ n and £2,7 are the energies of the photons 
in the Lab frame. The invariant diphoton mass must lie 
in the range 500 < m 77 < 575MeV/c 2 . For 773^, the pion 
candidates are fit to a common vertex to form an 77 can- 
didate. In the fit the 77 candidate mass is constrained to 
the nominal 77 mass, while the invariant three-pion mass 
is required to lie in the range 535 < < 565 McV/c 2 . 

B candidates are reconstructed from a hadron candi- 
date (tt + , 7T°, 77 77 , or 773^) and a dilepton candidate (e + e~ 
or p + pr). The hadron and dilepton candidates are fit to 
a common vertex, and the entire decay chain is refit. We 
make use of two kinematic, Lorentz-invariant quantities, 



tties and AE, defined as 

m ES = ^( S /2+p B -p ) 2 /E 2 -p% (1) 
AE = {2q B q Q -s)/2yTs (2) 

where y/s = 2E? is the total energy of the e + e~ sys- 
tem in the center of mass (CM) frame, q B and = 
(Eq,pq) are the four- vectors representing the momentum 
of the B candidate and of the e + e~ system, respectively, 
and ps is the three-momentum of the B candidate. In 
the CM frame, these expressions simplify to 

™ES = V^cam-PB 2 (3) 

AE = E* B - E* cam (4) 

where the asterisk indicates evaluation in the CM frame. 
These variables make use of precisely measured beam 
quantities. All B candidates arc required to have tties > 
5.1GeV/c 2 and -300 < AE < 250 MeV. The distribu- 
tions of these two variables are later fit to extract the 
ir£ + £~ and r)l + l~ branching fractions. 

A large background is present from B J/ipX and 
B —> tp(2S)X decays where J/ip and ip(2S) decay to 
l + (r . Here X represents a hadronic state, typically tt, 
rj, p, or K^*> . These events are removed from our data 
sample by rejecting any event with a value of mu consis- 
tent with originating from a J/tjj or ip(2S) decay. The 
rejected J/ip events are useful as they provide a con- 
trol sample which can be used to test the fit model. 
We also use these samples to estimate systematic uncer- 
tainties and to correct for differences between data and 
MC selection efficiencies. For the electron modes we re- 
ject events in the following regions about the J/ip mass: 
2.90 < mu < 3.20GeV/c 2 , or m ee < 2.90GeV/c 2 and 
AE < m ee c 2 — 2.875 GcV. For the muon modes the region 
is 3.00 < m w < 3.20GeV/c 2 , or m w < 3.00GcV/c 2 and 
AE < LHm^c 2 — 3.31 GeV. The rejection region about 
the ip(2S) mass is the same for electrons and muons: 
3.60 < m u < 3.75GeV/c 2 , or m u < 3.60GeV/c 2 and 
AE < muc 2 — 3.525 GeV. Introducing AE dependence 
into the region boundaries allows us to account for some 
of the effects of track mismeasurement and energy lost 
due to bremsstrahlung. 

The largest source of background comes from ran- 
dom combinations of particles from continuum events or 
semilcptonic B and D decays in BB events. Contin- 
uum events tend to be jet-like as the qq pair is produced 
back-to-back with relatively large momentum in the CM 
frame. The topology of BB decays is more isotropic as 
the B mesons are produced nearly at rest in the T(4S) 
rest frame. Semileptonic decays are characterized by the 
presence of a neutrino, e.g., missing energy in the event 
and non-zero total transverse momentum of the event. 
Due to the differences in these two background types we 
train separate artificial neutral networks (NNs) to reject 
each of them. By selecting inputs to the NNs which are 
independent of the final state we are able to train only 
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one NN for each lepton flavor. We do not train sepa- 
rate NNs for 7r + , 7T°, and r\. This increases the size of 
the training samples, improving the performance of the 
NNs. We train four NNs: one to reject BB background 
in thee + e~ modes, one to reject BB background in the 
p + p~ modes, one to reject continuum background in the 
e + e~ modes, and one to reject continuum background in 
the p + p~ modes. 

The signal training samples were assembled from equal 
portions of correctly reconstructed B + — > ir + £ + £~ , B° —> 
■k°£ + £-, B+ -> p+'£ + £~, B° -> p°£ + e~, B° -> rj^l+t-, 
and B° — > loI + £~ MC events. The size of the training 
samples, particularly the background training samples, 
that could be formed from events reconstructed as one 
of our signal modes was a limiting factor in the per- 
formance of the NNs. To increase the available statis- 
tics the other events from the p£ + £~ and lu£ + £~ modes 
were added to the training samples. The p + (p ) was 
reconstructed as 7r + 7r° (ir + ir~) with two-pion invariant 
mass m nn in the range 0.455 < < 1.095 GcV/c 2 

(0.475 < m„ n < 1.075 GeV/c 2 ). The uj was reconstructed 
as 7T + 7r~7r and required to have a three-pion invaraint 
mass lying within 50MeV/c 2 of the nominal uj mass. No 
r]3 V £ + £~ events were used in the training due to very low 
statistics for the background BB and continuum samples 
for these modes. For background we combined the MC 
data sets, either BB or continuum depending upon the 
classifier to be trained, from the six modes listed above 
and randomly select events from this data set to form the 
training sample. The performances of the NNs trained 
with samples from several different b — > dl + (r (global 
NNs) modes were compared with NNs trained specifi- 
cally for each of our four B — > tt£ + £~ modes (single mode 
NNs). The background rejection of the global NNs at a 
fixed signal efficiency was found to be similar to that of 
the single mode NNs. 

The input variables to the continuum NNs are related 
mostly to event topology and include the ratios of Fox 
Wolfram moments 
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the cosine of the polar angle of 
the thrust axis [32| of the event; the cosine of the polar 
angle of the thrust axis of the rest-of-the-event (ROE), 
which consists of all particles in the event not associated 
with the signal B candidate; the momentum weighted 
polynomials L\ [33| computed using tracks and EMC 
clusters; the cosine of the polar angle of the B candi- 
date momentum; and the x 2 probability of the B can- 
didate vertex fit. The input variables to the BB NNs 
reflect the effort to reject background from semileptonic 
B and D decays and include toes and AE constructed 
from the ROE; total momentum of the event transverse 
to the beam; missing energy in the event; momentum of 
the ROE transverse to the beam direction; momentum of 
the ROE transverse to the thrust axis of the event; cosine 
of the polar angle of the B candidate momentum; and \ 2 
probability of the B candidate and dilepton candidate 
vertex fits. The NN outputs show only weak correlation 
with the fit variables toes an d AE. 

Figure 2, as representative of the several neural net- 
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FIG. 2: Outputs of (a) the e + e~ BB neural network for a 
sample of B + — > ir + e + e~ signal (solid red) and BB back- 
ground (dashed blue) MC events, and (b) the p + p~ con- 
tinuum neural network for a sample of B° — s- iv p + p~ sig- 
nal (solid red) and continuum background (dashed blue) MC 
events (color available online). For both (a) and (b) the signal 
and background distributions are normalized to equal areas. 



works, shows the output of the e + e~ BB NN for a sam- 
ple of signal and BB background 7r + e + e~ MC events. 
Also shown is the output of the p + p~ continuum NN for 
a sample of signal and continuum background ir° p + pT 
MC events. Requirements on the NN outputs arc opti- 
mized for each of our eight modes to produce the lowest 
branching fraction upper limit. A description of the op- 
timization procedure is given in Section HO 

Due to their similarity to signal, B — > K^*H + £~ decays 
constitute a background that mimics signal by peaking in 
either one or both toes and AE. The b —> s£ + £~ transi- 
tion occurs at a rate approximately 25 times greater than 
the SM b — > d£ + £~ rate, and due to particle misiden- 
tification and event misreconstruction, its contribution 
is expected to be of the same order as the n£ + £~ sig- 
nal in the BABAR data sample. In the charged pion 
modes, B + K + £ + £~ peaks in toes as ir + £ + £~ signal 
but in AE near — 70MeV due to the misidentification 
of the kaon as a pion. There arc also contributions from 
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B° —> K®£ + £~ , where one of the pions from the decay 
is missed, and from B — > K*(— > K + ir)£ + £~ , where the 
pion from the K* decay is missed. In the case of Kg£ + £~ , 
the remaining pion and the two leptons are reconstructed 
as 7r+£+£-. For K+tt)£+£-, the K+ is misiden- 

tified as 7r + and reconstructed with the two leptons as 
ir + £ + £~ . In both cases, the decays peak in jhes like sig- 
nal but at AE < — 140 McV due to the missing pion. For 
K*£ + £~ the AE peak occurs at even lower values due to 
the kaon misidentification. For B° — > ir°£ + £~ , there is a 
similar background from B° — > K%(— > ir°ir°)£ + £~ decays 
where one ir is reconstructed along with the lepton pair 
as 7r °£ + £~. These events produce a peak in toes in the 
same location as tt°£ + £~ signal, but peak at smaller val- 
ues of AE due to the missing pion from the decay. In all 



three cases (K + £ + £~ in tt+£+£~, K° 
K*(-> K + 7r)£+t 



-)£+£- and 



ir + £+£~, and K°(-+ ^^)£ + £- in 
") we include a separate component in the fit model 
to account for the corresponding contribution. 



Hadronic decays such as B + — > 7r + 7r~7r + , where two 
pions are misidentified as muons, peak in both toes and 
AE similarly to signal due to the relatively small differ- 
ence between the pion and muon masses. This hadronic 
peaking background is modeled by a component in the 
fit. A dedicated data control sample is used to deter- 
mine its normalization and shape. This sample is con- 
structed from events where one lepton candidate passes 
the muon identification requirements but the other does 
not. The events in these samples are weighted with par- 
ticle misidentification probabilities determined from con- 
trol samples in BABAR, data. Studies of MC samples in- 
dicate that this background is consequential only for the 
7r/i + /i~ modes. 

After applying all selection criteria there are sometimes 
multiple candidates within a given mode remaining in an 
event. This occurs for approximately 20-25% (35-40%) 
of ir + e + e~ and ir°e + e~ (r/ 77 e + e~ and r/^e + e^) candi- 
dates, and 5-10% (25-30%) of n+fi+^r and ttV~V 



For Tre+e" and r^e+e", there is an additional back- (v-yj^fJ- and r? 3w /i+/i ) candidates. There tend to be 



ground that originates from two-photon events, given by 
the process e + e~ e + e~jj (e + e~)qq where q is a u, 
, or s quark. The background is characterized by a small 
transverse momentum of the pion and a large lepton- 
lcpton opening angle 8g+£-. There is also a correlation 
between the polar angles of the electron, 9 e - , and of the 
positron, 9 e + . The e~ tends to be in the forward direction 
while the e + tends to be in the backward direction, con- 
sistent with the e + e~ beam particles scattering into the 
detector. Events of this type are rejected using the fol- 



more events containing multiple candidates in the e + e~ 
modes due to the brcmsstrahlung recovery. For instance, 
there may be multiple candidates arising from the same 
7r + e + e~ combination where the bremsstrahlung photons 
associated with the e + or e~ are different. 

To choose the best candidate we construct a ratio Cr 
from the BB and continuum NN classifier output distri- 
butions of the signal and background samples. The ratio 
Cr is defined as 



lowing requirements. For n + e + e 



Tr°e+e- 



and 77 77 e + e 



we require pjj ad > 750 MeV/ c and A tr k > 4 where pjj ad 
is the hadron momentum in the CM frame and A tr k is 
the number of charged tracks in the event. Additionally, 
for 7r + e + e~ we require £i, n eut < 1-75 GeV, coaOg+g- > 
-0.95, and 6 e - > (0.57 9 e + -0.7 rad) where E l<neut is the 
energy of the highest energy neutral cluster in the event 
in the Lab frame. Similarly, 7r°e + e~ candidates must 
satisfy e - > (0.64 6 e + — 0.8 rad), and ?7 77 e + e - candi- 
dates are required to have 9 e - > (0.6 6 e + — 0.55 rad) and 
cos# f +^- > —0.95. These criteria were determined by 
maximizing the quantity e/y/NsB, where e is the signal 
efficiency and Age is the number of events lying in the 
sideband region 5.225 < toes < 5.26GcV/c 2 in data. We 
assume that the two-photon background in the toes side- 
band occurs similarly to the two-photon background in 
the region toes > 5.26GeV/c 2 . The optimization was 
carried out with all other selection criteria applied, in- 
cluding those on the NN outputs. 

To guard against possible background from B — > Dir 
and B — > Drj decays where D — > Kir, 7T7t, or rjir and the 
kaon or pions are misidentified as leptons, we assign the 
lepton candidates cither a kaon or pion mass and discard 
any event with a combination of /j^tt, or fj^r] 

with invariant mass in the range (1.83-1.89) GeV/c 2 . The 
probability of misidentifying a hadron as an electron is 
negligible, and this requirement is therefore only applied 
to the fi + fi~ modes. 



{V7eW + ^(y)) + V>Z§(x) + 

(5) 

where 'P S g S g( x ) C^cont^)) is the probability that a sig- 
nal candidate has a BB (continuum) NN output value 
of x (y). The quantities V g g(x) and V co ^ t (y) are de- 
fined analogously for background events. Signal-like can- 
didates have values of Cr near 1 while more background- 
like candidates have values near 0. If multiple candidates 
are present in an event, we choose the candidate with the 
greatest value of Cr as the best candidate. For events 
containing multiple candidates, this procedure chooses 
the correct candidate approximately 90-95% of the time 
for tt£+£- and 75-80% of the time for rj£+£-. The ratio 
Cr is used only to select a best candidate. 



IV. BRANCHING FRACTION MEASUREMENT 
AND UPPER LIMIT CALCULATION 

Branching fractions are extracted through an un- 
binned extended maximum likelihood fit to toes and 
AE with the fit region defined as toes > 5.225 GeV/c 2 
and -300 < AE < 250 MeV. The probability density 
functions (PDFs) in the fit model contain several com- 
ponents corresponding to the different contributions in 
the data set. To model the various components, we use a 
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combination of products of onc-dimcnsional parametric 
PDFs, two-dimensional histograms, and two-dimensional 
non-parametric shapes determined by a Gaussian ker- 
nel density estimation algorithm (KEYS PDF) [H[ • For 
components that are described by the product of one- 
dimensional PDFs, we are allowed to use such a model 
because Toes and AE are uncorrelatcd for these compo- 
nents. 



A. B+ n + £ + £- 

The ir + £ + £~ fit model involves four components: sig- 
nal, K + £ + t~ background, K° s /K*£ + £~ background, and 
combinatoric background. There is an additional compo- 
nent in B + — > 7r + fi + representing the B + — > 7r + 7r + 7r~ 
hadronic peaking background. The K + £ + £~ background 
arises from decays where the kaon is misidentified as a 
pion. The K + misidentification rate is such that the 
]{_+£+£- background in ir + £ + £~ is approximately the 
same size as the expected SM ir + £ + £r signal. Since the 
K + misidentification probability is well measured, it is 
possible to measure this background contribution directly 
from our data. This is done by simultaneously fitting two 
data samples, comprised by the B + —> tt + £ + £~ candi- 
dates and the B + — > K + £ + £~ candidates in our data 
set. The K + misidentification background to B + — > 
K + £ + £~ is included in the fit at a level fixed to the B + 
— > K + £ + £~ yield using the known misidentification prob- 
ability (which depends on the momentum of the kaon). 
The B + K + £ + £~ branching fraction that is measured 



from the simultaneous fit of the B 



+ 



-£+£- and B+ 



—> K + £ + £~ data samples provides an additional valida- 
tion of our procedure, since this branching fraction has 
been previously measured [37j . 

The K + £ + £~ sample is selected in exactly the same 
way as the k + £ + £~ sample except the charged pion iden- 
tification requirements are reversed and the J/tp and 
ijj(2S) rejection window includes the following regions: 
m ee > 3.20GeV/c 2 and l.llm ee c 2 - 3.67 < AE < 

surrounding the J/ip mass, 



m ee c 2 - 2.875 GcVfor n+e 



+ r, + . 



> 3.20GcV/c 2 and l.llm MM c 2 - 3.614 < AE < 
m^^c 2 — 2.925 GeV for w + fj, + p,~ surrounding the J/ip 
mass, and m« > 3.75GeV/c 2 and 1.11to«c 2 — 4.305 < 
AE < muc 2 — 3.525 GeV for both modes surrounding 
the -0(25") mass. Also, the AE window is -200 < AE < 
250 McV for K+e+e~ and -100 < AE < 250 McV for 

The tt+£+£- and K+£+£~ background m E s and AE 
distributions arc modeled by products of one-dimensional 
PDFs. The tt + £ + £~ signal and K + £+£~ background m ES 
distributions are described by a Crystal Ball function 
[35} . The 7r + e + e~ AE signal distribution is modeled 
by the sum of a Crystal Ball function and a Gaussian 
which share a common mean, while the 7r + /j + /i~ signal 
and both the K + e + e~a,nd K + fi + fj,~ AE distributions are 
modeled by a modified Gaussian with tail parameters 



whose functional form is given by 

(AE 



f(AE) = exp 



2<JL,ROt L ,R + «L,fl(A£ - fi) 



(0) 



where ctl and (ctr and an) are the width and tail pa- 
rameters used when AE < fi (AE > /1), respectively. A 
two-dimensional histogram models the contribution from 
B —> K2./K*£ + £~ decays. Combinatoric background is 
described by the product of an ARGUS function [3(| in 
toes with endpoint fixed to 5.29GeV/c 2 and a second- 
order polynomial in AE. The 7r + /x + /i~ hadronic peaking 
background component is modeled by a two-dimensional 
KEYS PDF H3. 

The PDF fit to the K + £ + £~ sample contains a simi- 
lar set of components. Signal K + £ + £~ distributions are 
modeled by the product of a Crystal Ball function in 
toes and the line shape of Eq. 6 in AE. The contri- 
bution from other b — > s£ + £~ decays is dominated by 
B K*(K + n)£ + £~ where the pion is lost. We use 
a two-dimensional histogram to model this background. 
Combinatoric background is modeled by the product of 
an ARGUS distribution in toeSj and by an exponential 
function for i^ + e + e~and a second-order polynomial for 
K+fi+fx-in AE. A KEYS PDF models the hadronic 
peaking background in K + /i + /i~ . 

In both the ir+£+£- and K + £ + t~ PDFs, the sig- 
nal and combinatoric background yields float along with 
the shapes of the combinatoric background PDFs. The 
K + £ + £~ background yield in the tt + £ + £~ sample is con- 
strained so that the B + — >• K + £ + £~ branching fractions 
measured in the ir + £ + £~ and K + £ + £~ samples are equal. 
All fixed shapes and yields are determined from exclu- 
sive MC samples except for the hadronic peaking back- 
ground which uses a data control sample. Normalizations 
of the K° s /K*£+£- component of the n+£+£- PDF and 
of K*£+£~ component in the K+£+£~ PDF arc fixed 
from efficiencies determined from MC samples and world 
average branching fractions (37j . 



B. B° -> tt°£+£- 

The B° — > tt°£ + £~ signal distribution is modeled by 
the product of a Crystal Ball function in toes and by 
the line shape given in Eq. [S] in AE. Background 
from B° -> K° s (-^ n Tr°)£ + £- decays is modeled by a 
two-dimensional histogram. The product of an ARGUS 
shape in Toes with an exponential function in AE mod- 
els the combinatoric background distribution. As in the 
7r + /i + /i~ and K + n + n~ PDFs, there is an additional com- 
ponent in the ir°ii + fi~ fit model devoted to hadronic 
peaking background which is described by a KEYS PDF. 

In the fit, only the signal ir°£ + £~ and combinatoric 
background yields along with the ARGUS slope parame- 
ter and argument of the exponential float. The signal and 
K®(— > Tr TT°)£ + £~ shapes are determined from fits to MC 
samples, and the K°(—^ ir°ir°)£ + £~ normalization comes 
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from efficiencies taken from MC samples and world aver- 
age branching fractions [37j . The shape and normaliza- 
tion of the peaking hadronic component are determined 
from a data control sample. 

C. B° -> rji+£- 

The rj£ + £~ fit model is simple, consisting of only three 
components, and is the same for all four rj£ + £~ chan- 
nels. The signal component is modeled by the product 
of a Crystal Ball function in tties & n d the line shape of 
Eq. 4 in AE. We include a component for events con- 
taining a signal decay where the signal B is incorrectly 
reconstructed, which we refer to as self-cross-feed. In 
these events the signal decay is typically reconstructed 
as a combination of particles from the B decaying to our 
signal mode and the other B. In most self-cross-feed 
events the dilepton pair is correctly reconstructed and 
the hadron is misreconstructed. The self-cross-feed con- 
tribution is represented by a two-dimensional histogram 
and its normalization is a fixed fraction of the signal 
yield with the fraction determined from signal MC. The 
self-cross-feed-to-signal ratio varies from 0.1-0.15 for the 
?y 77 channels to 0.25-0.3 for the rj^ channels. Combi- 
natoric background is described by the product of an 
ARGUS function in 7ties & n d an exponential function in 
AE. From studies of MC samples, we find no indica- 
tion of potential peaking background contributions from 
b s£ + £~ decays or any other sources. The i] 71 £ + £~ 
yield and the r\ i i- K £ + £~ yield are constrained in the fit to 
be consistent with the same B° — > r)£ + £~ branching frac- 
tion. The signal yield, combinatoric background yield, 
ARGUS slope and exponential argument float in the fit. 
All other parameters are fixed from MC samples. 



D. Lepton-flavor averaged and isospin averaged fits 

In addition to branching fraction measurements and 
upper limits for the B — > ir£ + £~ and B° — > r]£ + £~ 
modes we also present lepton-flavor averaged, isospin av- 
eraged, and lepton-flavor and isospin averaged results. 
The lepton-flavor averaged measurement of B(B + — > 
tt + £ + £~) is the branching fraction obtained from a si- 
multaneous fit to the 7r + e + e~ and 7r + /x + /x _ samples sub- 
ject to the constraint B(B + -> 7r+e + e~) = B(B + -> 
7r + /i + /i _ ). Here we have neglected the difference be- 
tween the electron and muon masses. The measure- 
ments of B(B° tt°£ + £-) and B{B° -4 rjl+i') are 
subject to a similar set of constraints and are deter- 
mined in an analogous way. The isospin averaged branch- 
ing fraction B(B — > Tte + e~) is the measured value 
of B(B + — > 7r + e + e~) after simultaneously fitting the 
7r + e + e~ and 7r°e + e~ samples subject to the constraint 
B{B+ Ti+e+er) = (t b o/2t b+ )B(B° 7r°e + e -) 
where t b o and Tg+ are the mean lifetimes of the neutral 
and charged B mesons, respectively [37| . An analogous 



expression is applied for the B{B — > tt(j, + /i ) measure- 
ment. The lepton-flavor and isospin averaged measure- 
ment of B(B -> n£+£-) is the value of B(B+ -> n+i+t') 
determined from a simultaneous fit to all four samples 
subject to both the lepton flavor and isospin constraints 
listed above. 



E. Upper limit calculation 

We set upper limits on the branching fractions follow- 
ing a method which utilizes the profile likelihood. Upper 
limits at the a confidence level (CL) are set by scanning 
the profile likelihood A as a function of the signal branch- 
ing fraction to determine where —2 In A changes by a per- 
centile of a x 2 random variable with one degree of free- 
dom. For a = 0.9 we look for a change in —2 In A of 1 .642. 
If the measured branching fraction is negative, we begin 
our scan from zero rather than the minimum [38j . This 
is a conservative approach that always produces physi- 
cal, i.e., non- negative, upper limits, even in the case of 
low statistics. Systematic uncertainties are incorporated 
into the limit by convolving the profile likelihood with a 
Gaussian distribution whose width is equal to the total 
systematic uncertainty. 

V. OPTIMIZATION OF SELECTION 

We simultaneously optimize the selection criteria for 
the two NN outputs and the PID selection criteria for the 
charged pions and leptons. BABAR employs algorithms 
which use outputs from one or more multivariate clas- 
sifiers to identify charged particle species. A few (3-6) 
standard selections on the outputs of these algorithms 
are used to identify particles of a given species with dif- 
ferent efficiencies. Greater identification efficiencies typ- 
ically imply greater misidentification rates. Due to this 
trade-off, it is not clear a priori which selection is best for 
a particular analysis. Therefore for each charged particle 
type (e - , tt + ) we optimize the PID requirements for 
the leptons and pions along with the NN output criteria. 

For the optimization we assume that B — > tt£ + £~ and 
B° —> rj£ + £~ occur near the center of the branching frac- 
tion ranges expected in the SM. Under this assumption, 
no statistically significant signal is expected, and the se- 
lection is optimized to produce the smallest branching 
fraction upper limit. We divide the BB and continuum 
NN output space into a grid and generate 2,500 paramet- 
rically simulated data sets per grid point according to our 
fit model. Each simulated data set is fit, and a branch- 
ing fraction upper limit is calculated. The figure of merit 
(FOM) for each point is the average branching fraction 
upper limit determined from the 2,500 data sets, and we 
take the combination of PID and NN output selection 
producing the smallest FOM as optimal. 

The results of the optimization show that the upper 
limits are rather insensitive to the PID selection. Also, 
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in the two-dimensional NN output space, there is a re- 
gion about the optimal selection where the FOM changes 
slowly, giving confidence that our optimization procedure 
is robust because the expected limits do not depend crit- 
ically on the NN selection requirements. 

The e + e~ modes use the same electron (muon) 

selection, while more efficient charged pion selection is 
favored for 7r + e + e~ and 'q^ 7T e + e~ than ir + [i + p~ and 
V37r^ + H~ ■ Tighter selection is favored on the continuum 
NN output than the BB NN output. The optimization 
favors looser requirements for the r)£ + £~ modes as the 
size of the background in these channels is much smaller 
than for tt£ + £~. 



VI. FIT VALIDATION 

We validate our fit methodology in three ways: (1) 
generating an ensemble of data sets from our fit model 
and fitting them with the same model ("pure pseudo- 
experiments"), (2) generating and fitting an ensemble of 
data sets with signal events from the BABAR MC simu- 
lation embedded into the data set ("embedded pseudo- 
experiments"), (3) extracting B — > J/ijiir and B° — > J/tpr] 
branching fractions from the BABAR data sample. 

From our studies of both pure and embedded pseudo- 
experiments, we find no significant source of bias in our 
fit. Distributions of branching fractions and their errors 
obtained from fits to these data sets are consistent with 
expectations. 

Measuring the B+ J/tpw+, B° J/i/jtt , B q 
J/tpr], and B + — > J/ifiK + branching fractions in the con- 
trol sample of vetoed charmonium events allows us to val- 
idate our fit methodology on data. We employ the same 
fit model to extract these branching fractions as we do 
for the ir+£+£~ } n°£ + £-, K+£+£~, and r]i + £- branching 
fractions. Fixed shape parameters and yields are deter- 
mined through fits to exclusive MC samples. We find 
that all measurements are in good agreement with world 
averages {37^ 

VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are included in the branching 
fraction upper limit calculation by convolving the pro- 
file likelihood with a Gaussian whose width is equal to 
the total systematic uncertainty. The systematic uncer- 
tainties are divided into "multiplicative" uncertainties, 
which scale with the true value of the branching fraction, 
and "additive" uncertainties, which are added to the true 
value of the branching fraction, independent of its value. 

A. Multiplicative uncertainties 

We list the sources of multiplicative systematic uncer- 
tainty below and their assigned values for each of the 



ir£ + £~ and r\£ + £~ signal modes in Table I. 

The systematic uncertainty in the measured number of 
BB pairs is estimated to be 0.6% (39| . 

The difference between the ir° reconstruction efficiency 
in data and MC has been studied in t + t - decays where 
one t decays via the channel r ± — s- e^vv and the other 
r decays via — > n v or — > p^v with p^ recon- 
structed as it^tt . The — > p^u yields are roughly pro- 
portional to the product of the 7r and 7r° reconstruction 
efficiencies, while the — > tx^v yields are proportional 
to the 7r reconstruction efficiencies. A correction pro- 
portional to the ratio of the r — > pv to r — > ixv yields in 
the data and MC samples is applied to better reproduce 
the data reconstruction efficiency in MC simulation. The 
uncertainty due to this correction is estimated as 3.0% 
per 7T°. We take the uncertainty in the ?7 77 reconstruc- 
tion efficiency associated with this correction to also be 
3.0% per ■q 11 . 

A correction to the MC tracking efficiency was devel- 
oped from the study of t + t~ decays where one r has 
a single charged daughter (1-prong decays) allowing the 
event to be identified as a t + t~ event and the other r has 
three charged daughters (3-prong decays). By measuring 
the event yields where the 3-prong r has either two or 
three tracks reconstructed, the track reconstruction effi- 
ciency can be measured. This efficiency can be used to 
correct the MC to match the efficiency measured in data. 
The systematic uncertainty associated with this correc- 
tion is estimated to be 0.3% per charged track taken to 
be 100% correlated among tracks in the event. 

We correct for the difference between the lepton PID 
selection efficiencies in data and MC by measuring the 
B + J/tpK + yields in data and J/tpK+ MC control 
samples with and without the PID selection requirements 
applied to both leptons. The ratios of the yields are 
used to correct the lepton particle identification selection 
efficiency derived from MC to match data. The error 
on the correction is taken as the associated systematic 
uncertainty, which ranges from 1.3-1.5%. The available 
statistics in the samples used to calculate the correction 
determine the size of the error which is associated with 
it. 

In an analogous procedure, we correct for the differ- 
ence between the charged pion PID selection efficiency 
obtained by measuring signal yields in high statistics 
B° J/tpK*°(^ K-JT+) data and exclusive MC con- 
trol samples with and without pion PID selection criteria 
applied. A correction is derived and the error on the cor- 
rection is taken as the associated systematic uncertainty. 
These uncertainties are approximately 2.5% and 3.5% for 
e + e~ and p^ pT modes, respectively. 

The high statistics of the B + J/ipK+ data and MC 
control samples are again exploited to derive a correc- 
tion for the NN output selection efficiency on MC. The 
J/tpK + signal yields were measured with only the BB 
NN output selection applied, only the continuum NN out- 
put selection applied, and with both selections applied. 
The error on the correction is taken as the associated 
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TABLE I: Multiplicative systematic uncertainties for the ir£ + £~ and rj£ + £~ modes. The lepton and tt PID and NN output 
selection efficiency correction uncertainties are determined using J/tp K^*' control samples, while the tracking and 7r°/?7 77 
efficiency correction and B counting uncertainties are taken from dedicated BABAR, studies. The total uncertainty is the sum 
in quadrature of the individual uncertainties. 





7r + e + e 




7T + /! + /l 


TV )1 + jj, 






77 77 /i+/j" 






0.6% 


0.6% 


0.6% 


0.6% 


0.6% 


0.6% 


0.6% 


0.6% 


7r°/?7 77 eff. 




3.0% 




3.0% 


3.0% 


3.0% 


3.0% 


3.0% 


Tracking eff. 


0.9% 


0.6% 


0.9% 


0.6% 


0.6% 


1.2% 


0.6% 


1.2% 


lepton PID 


1.3% 


1.3% 


1.4% 


1.4% 


1.4% 


1.4% 


1.5% 


1.5% 


PID 


2.5% 




3.5% 






2.3% 




3.7% 


NN cut 


1.4% 


1.3% 


1.5% 


1.5% 


1.4% 


1.3% 


1.5% 


1.5% 


Wilson coeff. 


2.7% 


2.3% 


1.0% 


1.9% 


3.1% 


3.1% 


0.3% 


0.9% 


FF model 


9.1% 


7.7% 


0.7% 


7.1% 


3.4% 


1.3% 


0.2% 


1.6% 


Total 


10.1% 


8.8% 


4.4% 


8.2% 


5.9% 


5.6% 


3.8% 


5.7% 



systematic uncertainty and ranges from 1.3-1.5%. 

We conservatively vary the Wilson coefficients CV, Cg, 
and C10 from their nominal values of —0.313, 4.344, and 
—4.669, respectively, by a factor of ±2 (e.g., C7 is var- 
ied to —0.157 and —0.616) and generate new simulated 
samples with all possible combinations of the varied Wil- 
son coefficients. For each varied sample we apply the 
full event selection and calculate the efficiency for that 
set of Wilson coefficients, taking the largest relative dif- 
ference between the varied Wilson coefficient efficiencies 
and the efficiency of our default model as the associated 
systematic uncertainty. 

Simulated MC samples using several different form- 
factor models were generated. Ultimately, we compare 
the efficiency from our default model with the efficiency 
calculated from the "Set 2" and "Set 4" form-factor mod- 
els of Ref. [!(|. The maximum relative difference be- 
tween our default model efficiency and the efficiency ob- 
tained with the "Set 2" and "Set 4" form- factor models 
is taken as the associated systematic uncertainty. There 
is a large variation in this uncertainty from one mode 
to another. The source of this effect is due in part to 
the correlation between the NN output and q 2 . Selec- 
tion on the NN outputs is mode dependent and therefore 
changes the q 2 dependence of the efficiency. Also, for the 
modes 7r + e + e~, 7r°e + e~, and r]~ n e + e~ we require that 
the hadron momentum be greater than 750 MeV/c in the 
CM frame. The hadron momentum is highly correlated 
with q 2 . Removing events with small hadron momen- 
tum also removes events with large q 2 . Differences at 
large q 2 between the differential branching fractions cal- 
culated using the default and alternative models lead to 
greater sensitivity to the choice of form- factor model, and 
therefore larger uncertainties associated with the choice 
of model. 

The uncertainty in the efficiency due to the size of the 
simulated MC samples is less than 0.1% and is negligible. 



B. Additive uncertainties 

We consider the following sources of additive system- 
atic uncertainty with their values given in Table II. 

The fixed parameters of the ir£ + £~, r/£ + £~, and 
K + £ + £~ signal and the K + £ + £~ background PDFs are 
varied individually within the errors obtained from fits 
to exclusive MC samples, and the data sample is re- fit. 
For simultaneous fits we additionally vary the efficiencies 
within their uncertainties, and for rj£ + £~ we vary the 
self-crossfeed-to-signal ratio by ±10%. The size of the 
variation is arbitrary but conservative enough since the 
number of expected self-crossfeed events is at most 0.15. 
The difference between the branching fraction from this 
fit and that from the nominal fit is taken as the associ- 
ated systematic uncertainty. We take the largest change 
in the branching fraction as the systematic uncertainty 
associated with each fixed quantity. The uncertainties 
from individual variations are summed in quadrature. 

Non-parametric PDFs include the two-dimensional his- 
tograms and KEYS shapes. We vary the binning of the 
two-dimensional PDFs, increasing and decreasing them 
by a factor of two. The data sample is re-fit, and we take 
the largest change in the branching fraction as the asso- 
ciated systematic uncertainty. For the KEYS shapes, we 
increase and decrease the width of the Gaussian kernel 
used to generate the shapes and take the largest change 
in the branching fraction as the associated systematic un- 
certainty. If there are multiple non-parametric PDFs, we 
add their associated uncertainties in quadrature. 

The hadronic peaking background yields arc fixed from 
the control sample of hadronic decays and are varied 
within their statistical uncertainties. The data sample 
is re- fit, and we take the largest change in the branching 
fraction as the associated systematic uncertainty. 

We fix the K°/K*£ + £~ yield in the ir+£+£- PDF, the 
K*£+£~ yield in the K+£+£~ PDF, and the K%£ + £~ 
yield in the tt°£ + £~ PDF. These values are determined 
from efficiencies taken from exclusive MC samples and 
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TABLE II: Additive systematic uncertainties for the B — > ty£ + £ and B° — > r/£ + £ channels. The total uncertainty is the sum 
in quadrature of the individual uncertainties. All uncertainties are given in units of 10 -8 . 



Mode 


7r + e + e 


TT°e+e- 


7T + /! + /l 


o + . - 

7T fl fl 


r/e e 


rffi+fi 


Fixed parameters 


1.9 


0.2 


0.6 


0.3 


0.6 


0.4 


Non-parametric shapes 


< 0.1 


< 0.1 


0.7 


0.5 


0.1 


0.1 


Hadronic peaking bkg yields 






< 0.1 


< 0.1 






Non-hadronic peaking bkg yields 






< 0.1 


< 0.1 






Total 


1.9 


0.2 


0.9 


0.6 


0.6 


0.4 



the current world average branching fractions for these 
modes 37[ . We vary the yields according to the errors on 
their branching fractions and re-fit the data sample. The 
change in the branching fraction from its nominal value 
is taken as the associated systematic uncertainty. In Ta- 
ble II, these uncertainties are classified as "Non-hadronic 
peaking bkg yields" . 



VIII. RESULTS 

We extract branching fractions by fitting the data set 
with the fit model described in Section IIVI Projections 
of the PDFs and data sets in tties and A IS are shown for 
the isospin averaged B — > ire + e~ fit, the lepton-flavor 
averaged B° — > r)£ + £~ fit, and the isospin and lepton- 
flavor averaged B — > tt£ + £~ fit in Figs. 3-5, respectively. 
Figure 6 shows —2 In A as a function of the branching 
fraction for the tt£+£-, n + i + £~, tt°£+£- t]£+£-, rje+e^, 
and r7/x + /x _ measuremcnts. Branching fraction measure- 
ments and upper limits at 90% CL are given in Tabic III 
for each mode. 

As a cross-check, we measure the B + — > K + £ + £~ 
branching fractions and find them consistent with the 
current world averages j37| . 

We set upper limits on the lepton-flavor averaged 
branching fractions of 



B{B + -> Ti + £ + r) < 6.6 x icr 8 , 
B(B° -> ir°£+£-) < 5.3 x 10~ 8 , 
B(B° -> r]£ + £-) < 6.4 x 10~ 8 , 



(7) 
(8) 
(9) 



all at the 90% CL. A lepton-flavor and isospin averaged 
branching fraction upper limit of 



B(B -> ir£ + r) < 5.9 x 10" 



(10) 



is set at 90% CL. Branching fraction measurements and 
upper limits at 90% CL for the modes B+ 



B u -> 7r u e+e 



rje + e , and B 



B+ -> TT+fX+fJ,-, B° -> Tr°ti+fj,-, B" -> 
i]fi + fi~ arc listed in Tabic III. 
In conclusion, we have searched for the rare decays 
B -> n£ + £~ and B° -> rf£ + £~ in a sample of 471 million 
BB decays and observe no statistically significant signal 
in any of the decay channels studied. We set a lepton- 
flavor and isospin averaged upper limit at the 90% CL 



of B{B -> tt£+£-) < 5.9 x 10" 8 , within a factor of three 
of the SM expectation. We also set lepton-flavor aver- 
aged upper limits of B(B+ -> 7r + £ + £r) < 6.6 x 10" 8 and 
B(B° -> n°£ + £-) < 5.3 xlO" 8 . Branching fraction upper 
limits at 90% CL have also been calculated for the modes 
B+ -> Tr+e+e-, B° -> n°e + e-, B+ -> tt+ fi+ pT , and 
B° —> 7r°/i + /i _ . Our upper limits on the B° —> ir°e + e~, 
B° —> 7r°/i + /i _ , and B° — > ir°£ + £~ branching fractions 
are the lowest upper limits to date. The results pre- 
sented for the ir£ + £~ modes supersede those of the pre- 
vious BABAR analysis [2(J. We have also performed the 
first search for the decays B° —> r/£ + £~ and set an upper 
limit on the lepton-flavor averaged branching fraction of 
B(B° -> r)£+£~) < 6.4 x 10~ 8 at the 90% CL. Upper 
limits at 90% CL for the B° -> jje+e' and B° 
branching fractions have been reported. 
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TABLE III: B — > tt£ + £~ and B — > r)£ + E~ efficiencies, yields, branching fractions, and branching fraction upper limits at the 
90% CL. The error on the yield is statistical. The first error quoted on the branching fractions is statistical while the second 
is systematic. Branching fraction upper limits include systematic uncertainties. 

Mode e Yield g(10~ 8 ) Upper Limit (10~ 8 ) 

B+^-K+e+e' 0.199 ^-tl 4.3^ ? ± 2.0 12.5 

B°^iv e + e- 0.163 1-OIlI 1.2i^o ± 0.2 8.4 

B° -> r)e + e~ -4.0tg° ± 0.6 10.8 

B° r7 77 e+e- 0.164 — 1.2±|;| 

B° -> 773,, e+e" 0.115 — 0.5±^;g 

B+^-K+ft+Li- 0.140 — 0.5±|;| -0.6±£a±0.9 5.5 

B ^ttV + m" 0.115 -0-2±o;° -1.0±|;2±0.6 6.9 

B°^rui + fj,- -2.0lg;g±0.4 11.2 

S° -> ?? 77 /J + /i~ 0.102 -0.4±i;3 

S° -4 rtsn^n" 0.063 -O.lt^.I 

B -> 7re+e~ 4.0±!^ ± 1.6 11.0 

S -> 7r/j+/i~ -O.Qt^o ± L2 5 -° 

B+ -»■ tt+^+zT 2.5±|:| ± 1.2 6.6 

B°^tt°£ + £- l-2±|:i ± 0.2 5.3 

B° r]£ + £- -2.8±jj;f ± 0-3 6.4 

B -» ttI + 1~ 2.5±| g ± 1.0 5J) 
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B+ -> Tr+e+e" B+ ->■ Tf+e+e- S° -> TrV+e" 
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FIG. 3: Fit projections of m^s (top) and AE (bottom) for the isospin averaged B — ► 7re + e~ fit to the 7r + e + e _ ((a) and (b)), 
K + e + e~ ((c) and (d)), and ir°e + e~ ((e) and (f)) data sets. Points with error bars represent data. The curves are dotted 
(magenta) for combinatoric background, dot-dashed (gray) for K* /K%£ + &~ background, dashed (green) for K + £ + £~ signal 
and background, and solid (red) for tv£ + £~ signal. The solid blue curve represents the total fit function, (color available online) 
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FIG. 4: Fit projections of jtibs (top) and AE (bottom) for the lepton-flavor averaged fit to the r\ 11 e r e~ ((a) and (b)), ?73,re + e _ 
((c) and (d)), r; 77 /i + ^t~ ((e) and (f)), and ^3^/x + \iT ((g) and (h)) samples. Points with error bars represent data. The curves 
are dotted (magenta) for combinatoric background, dashed (green) for self-crossfeed, and solid (red) for r]£ + £~ signal. The 
solid blue curve represents the total fit function, (color available online) 
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B+ -> 7r + e + e" 



B+ -> K+e+e 
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FIG. 5: Fit projections of rriES and AE for the lepton flavor and isospin averaged ir£ + £ fit to the 7r + e + e ((a) and (d)), 
K+e + e~ ((b) and (e)), 7r°e + e~ ((c) and (f)), tt+ ^+ /x" ((g) and (j)), K + fi + fi~ ((h) and (k)), and 7r°^ + /i" ((i) and (1)) data 
sets. Points with error bars represent data. The curves are dotted (magenta) for combinatoric background, triple dot-dashed 
(cyan) for hadronic peaking background in the fi + modes, dot-dashed (gray) for K* /K®l + £~ background, dashed (green) 
for K + £ + £~ signal and background, and solid (red) for n£ + £~ signal. The solid blue curve represents the total fit function, 
(color available online) 
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FIG. 6: The negative logarithm of the profile likelihood as a function of branching fraction for (a) B — ¥ tv£ + £ , (b) B° — > r]£ + £ , 
(c) B + -5- n + £ + r, (d) B° 7?e+e", (e) B° % a £ + r, and (f) 5° -)• The solid, black lines are the negative log 

likelihood curves including only statistical errors, and the dashed, green lines correspond to the same curves convolved with a 
Gaussian distribution whose width is equal to the total systematic error, (color available online) 



